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ABSTRACT 

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric 
fiber composite, interdigitated electrode blade twist actuators is described. The analysis consists 
of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady 
aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial 
differential equations of motion to develop a system of ordinary differential equations suitable for 
dynamics simulation using numerical integration. The twist actuation responses for three 
conceptual full-scale blade designs with realistic constraints on blade mass are numerically 
evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant 
elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight 
conditions for interdigitated electrode poling configurations. Twist actuation for the 
interdigitated electrode blades is also compared with the twist actuation of a conventionally poled 
piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode 
actuators was found to be four to five times larger than that obtained with the conventionally 
poled actuators. 


1. INTRODUCTION 

Significant undesirable fixed system and rotating system vibratory loads continue to exist on all 
helicopters. These vibratory loads are primarily the result of responses generated by unsteady 
aerodynamic loads acting on the main rotor system while the helicopter is in forward flight. Since 
high vibratory loads result in increased maintenance requirements and poor ride quality, much 
effort has been devoted to finding means to eliminate or reduce them. 

Many mechanical devices, mounted either in the fixed system (fuselage), or rotating system (hub, 
and occasionally blades), have been developed for this purpose. 1, 2 The majority of these 
mechanisms essentially attempt to alleviate the undesirable vibratory loads by introducing 
counteracting inertial and damping forces. This approach has the drawback of requiring the 
addition of mechanisms, with attendant mass and complexity and sometimes performance 
degradation, to the vehicle. In contrast to these methods, the techniques of higher harmonic 
control (HHC), 3 4 and individual blade control (IBC), 5 ' 6 ’ 7 ’ 8 seek to eliminate vibrations at their 
aerodynamic source; by directly altering the unsteady aerodynamic forces acting on the rotor 
blades. This is accomplished by introducing cyclic variations in the root pitch of each blade, either 
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by clever manipulation of the existing (swashplate) control system, or through the use of auxiliary 
pitch actuators mounted within the rotating system. 

HHC active control through a swashplate-type mechanism is restricted by fundamental rotor 
dynamics to cyclic pitch inputs at the blade passage frequency (M2), and the next higher and 
lower harmonics of the blade passage frequency ((N ± 7)£2), only, where N is the number of rotor 
blades, and £2 is the rotational speed. HHC can also be limited by the typically low bandwidth 
characteristics of swashplate hydraulic actuation systems. 

The IBC approach, on the other hand, enables individual control of blade pitch at essentially 
arbitrary frequencies. This is accomplished through auxilliary pitch control actuators placed 
within the rotating system. With few exceptions, hydraulics based systems have proven to be the 
best means of providing the power and bandwidth in the rotating system necessary for useful IBC. 
However, the added complexity associated with delivering hydraulic power from the fixed system 
to the rotating system has made practical implementation of IBC on production helicopters 
prohibitive. 


1.1 Smart material IBC actuation schemes. 

Methods of using smart materials to deflect trailing edge control flaps, 9, 10, 11 or actuate blade 
twist ' electromechanically have been advanced as a more practical alternative means of 
implementing some form of IBC. Electrical power for smart material actuators has the advantage 
of being easily deliverable from the fixed system to the rotating system through conventional 
slipring devices. This avoids the complications associated with transferring fixed system hydraulic 
power for individual blade actuators to the rotating system. The bandwidth characteristics of 
certain smart material actuators are also superior to those of hydraulic and conventional 
electromechanical, or servomotor, actuators. 

Currently, trailing edge flap actuation has the advantage over other smart material IBC schemes in 
that available smart material power and displacement capabilities are very near those required for 
a practical rotor blade flap actuator design. Disadvantages of the trailing edge flap include the 
undesirable addition of complicated mechanisms to the rotor system, and to a lesser extent, the 
degradation of flap aerodynamic force and moment performance due to hinge lines. 

Blade twist actuation is attractive by virtue of being mechanically simple and aerodynamically 
efficient. However, as rotor blade structures tend to be relatively stiff in torsion, induction of 
useful levels of twist using smart materials can be difficult. In addition, the added mass of the 
embedded smart materials can also be prohibitive when developing practical induced twist blade 
systems. Despite these drawbacks, some encouraging developments in twist actuation of smart 
material structures continue to be made. Most recently, research in anisotropic twist actuation of 
plate structures using piezoelectric fiber composites (PFC) 14 ' 15 has demonstrated that relatively 
high levels of twist actuation are potentially achievable. The application of interdigitated 
electrode technology (IDE) 16 can in principle enhance the performance of these materials even 
further. 
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1.2 Previous work related to analysis of on-blade active control. 


To date, relatively little analytical work has been reported detailing the aeroelastic behavior of 
rotor blades incorporating either trailing edge control flaps or controllable twist. Millott and 
Friedman 17 have given the most comprehensive treatment so far with respect to the unsteady 
aerodynamics of rotor blade trailing edge control flaps, and their employment as vibradon 
reduction devices. Specific methods for actuation of control flaps were not addressed in their 
work however. Use of control flaps to generate higher harmonic controllable twist was 
investigated analytically and experimentally in full-scale wind tunnel tests performed in the 1970's. 
The analytical portions of this investigation were conducted by Lemnios and Dunn, 18 with wind 
tunnel results reported by Wei and Weisbrich. 19 Significant vibration reductions were reported 
using flap actuation at OQ, Id, 2d, 3d, and 4d frequencies. However, flap actuation was 
achieved only by use of a relatively complex system of mechanical linkages, which would not be 
generally suitable for flight test or production vehicles. A related investigation of the performance 
effects of variable camber was also carried out by Dadone, Cowan, and McHugh. 20 Analytical 
results suggested that significant improvements in overall rotor performance were achievable, 
especially when airfoil camber changes were used to induce blade twist. Techniques for affecting 
camber changes on rotor blades were discussed only in conceptual terms. 

Even less has been reported on the analysis of the aeroelastic behavior of rotor blade structures 
specifically incorporating embedded smart material actuators. Song and Lebrescu 21 developed the 
equations of motion for a rotating, thin-walled, cantilevered beam structure incorporating 
embedded piezoelectric actuators. No aerodynamic forcing was included in their study, and 
actuation of torsional motion was not considered. Nitzsche and Breitbach recently reported the 
results of a study to evaluate the ability of embedded piezoelectric materials to attenuate out-of- 
plane bending and torsional vibrations on a rotor blade structure. 22 To accomplish this, they 
developed a rotor blade aeroelastic model incorporating quasi-static aerodynamics and a 
"directionally attached piezoelectric crystal" bending-torsion actuation scheme similar to that 
developed by Barrett. They concluded that the lightly damped torsional blade modes could be 
significantly affected on a practical blade structure without saturation of the piezoelectric 
materials. 


1.3 Scope of this effort. 

For the most part, aeroelastic analysis of rotor blades incorporating on-blade actuation is still very 
much in its infancy. In particular, there is a lack of simple analytical models suitable for 
conducting preliminary conceptual control and design studies for embedded smart material blade 
structures. In light of this, and in order to gain greater insight into the control and aeroelastic 
response issues related to induced twist smart structure rotor blades, a simple aeroelasticity model 
for a piezoelectric fiber composite twist actuated helicopter rotor blade has been developed by the 
authors. This model is derived specifically for use in the investigation of phenomena related to 
torsional control and response of helicopter rotor blades incorporating piezoelectric twist. In this 
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paper, a description of the derivation and numerical implementation of this model is given. 
Additionally, numerical results demonstrating the twist actuation potential of three conceptual 
full-scale helicopter blade designs are shown. The work reported here is thought to be the first 
specifically related to the aeroelastic analysis of piezoelectric fiber composite twist actuated 
helicopter rotor blades. 


2. ANALYTICAL MODEL DESCRIPTION 

A rotor blade aeroelasticity model may be thought of as a unification of several basic analytical 
submodels. These submodels are, at a minimum, 1) a formulation of the blade structural 
dynamics, including all forces related to blade rotation, 2) a formulation of the aerodynamic forces 
and moments acting on the blade, 3) an inflow model, and 4) a control formulation. These 
submodels may be coupled in many ways depending on the degree of complexity of the overall 
formulation. In this development, they will be conceptually linked as shown in Figure 1. 
Theoretical descriptions of each submodel are given below. 


2.1 Structural formulation 

The equations of motion used here to describe the elastic torsion and out-of-plane bending 
behavior of an isolated helicopter rotor blade are adapted from the general elastic bending and 
torsion deformation equations developed by Kaza and Kvatemik. 23 Due to the complexity of 
these equations, and elastic rotor blade equations of motion in general, it is usually necessary to 
apply some simplifying assumptions to the complete set of equations in order to obtain a more 
mathematically manageable model. An ordering scheme approach is used here to accomplish this. 
Use of such a procedure enables one to avoid a great deal of the complicated algebra associated 
with retention of "high order" nonlinear terms in the equations, and ensures that the most 
physically significant terms are retained. 

The first step in employing an ordering scheme is to rewrite the equations of motion in a 
nondimensional form. An order is then assigned to each nondimensional term in the equations 
relative to an assumed scale factor, 6 , and the equations are rewritten retaining only terms of the 
lowest order, and terms of one order greater. Terms of two orders greater and higher are usually 
discarded. For this study e will be assigned a value of 0. 1 0, which is equal to the assumed order 
of magnitude of torsional deformations, hence, O(0)=l, O(e)= 0.I0, O(e 2 )=0.01, etc. 

Evaluation of the order of each nondimensional term is based on the intended application of the 
final equations, e.g., stability, vibration, etc., and representative full-scale helicopter structural and 
flight parameters. The ordering of parameters used in this study is based on schemes applicable to 
rotorcraft vibrations, and is shown in Table 1. 

Three additional assumptions related to the blade geometry are also made. These are, 1 ) that the 
blade precone angle and built in twist are assumed to be zero, 2) the blade structural cross-section 
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is assumed to be doubly symmetric, and 3) the blade pitch radius of gyration is assumed to be 
approximately equal to the k m2 cross section integral (i.e., k m i/k m2 «l). 

Applying the ordering scheme, with the additional assumptions given above, to the equations of 
Reference 23 yields the following partial differential equations of motion for blade out-of-plane 
bending and torsion. (Blade coordinates are illustrated in Figure Al.) 

Out-of-plane (flapwise) bending ( 0(e ), 0(£*)): 
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Torsional deformation (0(e 4 ,e 5 )): 
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The blade section tension force, T , is based on an 0(1 ) approximation, and is given by 

T = il 2 R 2 [ rrvcdx . 

Jx 

(3) 

Temporal differentiations (denoted by *) in the above equations are performed with respect to the 
rotor azimuth, v|/, which, for constant rotational velocity il, is \p = ilt . Spatial differentiations 
(denoted by + ) are conducted with respect to the nondimensional blade radial location, x = x/R . 
Descriptions of the coordinate systems used to describe the blade motions may be found in 
Appendix A. 

Note that the single underlined terms in Equation 1 are actually of 0(e 4 ), but are retained in order 
to ensure mathematical symmetry in the resulting mass and stiffness matrices. Double underlined 
terms in Equation 1 are also of 0(E 4 ), but are kept so that the equation is a more physically 
meaningful representation of out-of-plane bending. Technically, retention of all 0(e 4 ) terms 
would require inclusion of additional terms associated with the in-plane ("lead-lag") bending 
deformation of the blade; terms that have been consistently ignored in this treatment. A higher 
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order approximation of the tension force would be also be necessary for consistency. We re- 
emphasize that these equations were developed primarily for use in investigations related to blade 
torsional response and control. Consequently, studies of out-of-plane vibratory response using 
the present formulation should be done with some degree of caution. 

Qpe , in Equation 2, represents the piezoelectric induced twisting moment, which will be derived 
in the control formulation section below. 

A modified Galerkin procedure 24 is used here to obtain modal solutions to Equations 1-3. In this 
case, superposition solutions for w and (j> of the form, 

L 

/= l 

(4) 

♦(*.V) = XMv>M*). 

m=l 

(5) 


are assumed, where L and M are the number of out-of-plane bending and torsional modal 
functions respectively. In the modified Galerkin procedure, these modal functions need only 
satisfy the geometric boundary conditions on the blade. Work due to any nonfulfilled "natural" 
boundary conditions is accounted for with additional boundary terms in the equations. 
Substituting Equations 4 and 5 into Equation 1, and performing the appropriate integrations, 
yields a set of L ordinary differential equations of the following form: 




t*. f —~ y n , W"W„"dx 
tt Jo mQ 2 R i ' " 


+X W ' 

/=1 


1 T fi 

w; w: dx + f — w; + w n ++ dt + 

Jq r\2 n 4 ‘ n 


Jo mfrR 


2 WV I T n 


m£l 2 R A 

+ £ < ) > -" X §mf o ex<S> n W;dx 

m=l m —\ 

= f — ; w n dx - fe \'eW dx + 0 f exW+dx ) 

Jo mQ 2 R " l J o n Jo » ) 


Q, 2 R 2 I'mdx 

Jo 


W, + (0)W„ + (0) ’ 

) 


( 6 ) 

where n = 1, L . A similar procedure performed on Equation 2 yields an additional set of M 
ordinary differential equations; 
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The AT p and AT e terms in Equations 6 and 7, which do not appear in Equations 3 and 4, are used 

to represent finite stiffnesses present at x = 0 . These terms can be used to account for the 
stiffness of a mechanical flapping spring placed at the blade root, or the inherent flexibility of the 
pitch control system. 

It should be noted that all stiffness, damping, and inertia terms contained in the above equations 
represent the structural properties of the combined piezoelectric/passive material blade structure. 
The derivation of the stiffness terms is contained in Appendix B. 


2.2 Aerodynamic formulation 

The sectional lifting forces and moments are calculated using a technique based on the ONERA 
dynamic stall model developed by Tran and Petot. 25 The ONERA model uses differential 
equations in time to describe the unsteady aerodynamic lifting forces and pitching moments, 
including dynamic stall effects, acting upon an airfoil section undergoing arbitrary pitch and 
plunge motion. 


2.2.1 Section lift formulation 

Modifications to the ONERA model for general use in rotorcraft aerodynamic formulations have 
been made by Peters, 26 with nondimensional circulations used as the state variables instead of 
aerodynamic coefficients. The simplified lift circulation equations of Reference 26, which are well 
behaved in the reverse flow region of the rotor disk, but do not give lift reversal, are used here (8, 
9). 
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L 0 in equation (8c) is the nondimensional apparent mass lift. L x and L y are the components of 

the nondimensional lift in the airfoil section X and Y directions respectively. The nondimensional 

* 

section velocities U, U Y , and Uz, and section angle of attack, a, are derived in Appendix C. e in 
equations (9) is defined as the geometric rate of rotation of the airfoil with respect to the air mass 
(see Appendix C), and is given here by 


* de 
e s — 
d \ |/ 


= e con + ( t) +M, + - 


( 10 ) 


AC, in Equations 9 is the difference between the linear model static lift coefficient, C zt and the 
measured stalled lift coefficient, C a , i.e.. 


AC Z =C zt —C zs . 

Plots of the static lift coefficient data used in this model are shown in Figure 2. These curves are 
extrapolated from curves given in References 28 and 29. 

The angle of attack dependent coefficients (.v, A, 8, d, w, and e) are derived from experimental 
two-dimensional unsteady airfoil tests using the parameter identification scheme described in 
Reference 25. Parameter values for the ONERA OA212 rotorcraft airfoil are used in this model, 
and are shown in Table 2. 


2.2.2 Section Pitching Moment Formulation 
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Improvements to the basic ONERA pitching moment formulation have been made by Petot, 27 
with further modifications made by Chouchane 28 (given also by Peters, Chouchane, and Fulton* 9 ), 
and this is the representation used in this model. In this approach, the unstalled component of C m 
is given explicitly through the static moment coefficient, which is a function of angle-of-attack 
only. This results in the elimination of one state per spanwise aerodynamic evaluation point in the 
model. Plots of the static pitching moment data used here are shown in Figure 3. These curves 
are extrapolated from curves given in References 28 and 29, and from data provided by Tang. 

The stalled contribution to the section pitching moment is calculated using a circulation based 
model similar to that developed for section lift. The second order differential equation describing 
this stalled pitching moment circulation, defined as T m2 = UC m2 , is shown below (1 1). 


k ; r m2i + Cl_ k, r n 


'2i+r m ,r m2i 


= —rU i AC m . 

m i i m i 


E m ,k, U yi 


( 11 ) 


The coefficients a m , r m , and E m in Equation 1 1 were found by Petot (Reference 27) to have similar 
characteristics for many airfoils. Expressions for these coefficients, omitting subscripts, may be 
written as 


a m = a o +a 2 AC z 2 , 
r m =(r 0 +r 2 AC/) 2 , 

E m =E 2 AC 2 . 


(12a) 

(12b) 

(12c) 


Values of a 0 , a 2 , r 0 , r 2 , and E 2 used in the present formulation are taken from the generic "mean 
airfoil” values of Reference 27, and are provided in Table 3. 


2.2.3 Airloads calculation 

The aerodynamic forcing integrals present in Equations 6 and 7 are calculated by first evaluating 
the sectional aerodynamic forces and moments per unit length at N discrete points along the blade 
span. Section aerodynamic forces and moments are further assumed to be constant over the 
width of each section. Consistent with the ordering scheme assumed in the structural formulation 
above, the final expressions for these aerodynamic loading integrals are 
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r , in Equations 13 and 14, is the radial location of the inboard edge of the ith aerodynamic 
section. 


r i= x i ~ A i/2, 

where x t is the nondimensional radial location of the ith aerodynamic evaluation point, and A ( 
the associated nondimensional section width. 


2.3 Inflow model 


A uniform rotor inflow model, based on momentum theory (e.g., Gessow and Myers 31 ), is used in 
this formulation. A uniform inflow assumption, i.e., the assumption of constant inflow velocity at 
every location across the rotor disk, is not unreasonable for hovering or vertical flight conditions, 
although it is not very realistic for forward flight. Nevertheless, it is used here for computational 
simplicity in the numerical model. More complex inflow representations will be incorporated into 
future versions of this analysis. 


For an implicitly trimmed rotor operating at thrust coefficient C T , and forward flight advance ratio 
p, the uniform induced velocity ratio with respect to the tip-path-plane may be written as 


^tpp = ptana, + k i , 


where 


(15) 


X, = 


2-^)1 2 + X - 


TPP 


(16) 


An Newton-Raphson iterative approach, from Johnson 32 is used to numerically evaluate Xtpp in 
the above expression, i.e. (17) 
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Usually only two or three iterations are required for the calculation to converge, starting with 

C T 


Xjpp = |i tan a t + 


2> 2 + C r /2 ' 


(18) 


2.4 Control formulation 


2.4.1 Trim control 

Swashplate control angles are calculated using trim equations adapted from the harmonic balance 
equations of Reference 32. These equations, rewritten assuming zero blade twist and no first 
harmonic flapping with respect to the rotor shaft, are given by 
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where |5 0 in equation 19c is the estimated blade coning angle given by 
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The control settings predicted using these trim equations are reasonably effective at eliminating 
first harmonic flapping in unstalled flight. However, under heavily stalled forward flight 
conditions, these equations will fail to trim the rotor effectively. 


2.4.2 Piezoelectric twisting moment 

The piezoelectric actuator equation development detailed here borrows heavily from the approach 
used by Bent, Hagood, and Rogers. 15 Their equations, however, have been rewritten here 
explicitly for the case of in-plane polarization of the piezoelectric lamina material, or interdigitated 
electrode (IDE) poling. 16 Such an IDE configuration, in principal, allows the so-called "d 33 
effect" to be used to maximum advantage to enhance the inplane strain actuation capabilities of 
the piezoelectric lamina. This in turn should result in an increase in torsional actuation capability 
of the lamina with respect to conventionally poled schemes. 


Development of the actuator equations begins with an alternate version of the standard linear 
piezoelectric constitutive relations. 33 Assuming structural orthotropy in the piezoelectric material, 
and applying the plane stress assumption (T 3 =T 4 =T 5 =0) to the constitutive relations yields 


d; 


\t T 

c n 

0 

0 

d 33 

d 2 


0 

f t 

C . 22 

0 

0 

.a 


0 

0 

p r 

c 33 

0 



d 33 

0 

0 

s u 

*2 


*3. 

0 

0 

4 

A. 


0 

d u 

0 

0 


0 

0 


J 12 


'’22 

0 



The S 3 , S 4 , and S 5 strains, although not necessarily zero, will be neglected here. 


( 20 ) 


The poling direction of the piezoelectric material will be defined to be in the material 1-axis 
direction (see Figure 4), rather than in the more familiar out-of-plane direction (3-axis). The usual 
indexing nomenclature in the piezoelectric free-strain coupling coefficients will be retained here 
however. This is solely to allow the d 33 , and d /3 , piezoelectric material coefficients to be more 
readily identified in the resulting equations. (Note that the conventionally poled case may easily 
be obtained from these equations by setting d 33 equal to d 33 .) 
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Rewriting Equation 20 with strains (5) as independent variables: 


D, 


e n 
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0 

D, 


0 
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0 0 e 33 

e s 22 0 0 
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0 0 c,j 
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0 

0 


12 

c £ 

l 22 
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0 

e i5 

0 

0 

0 

4 


J 6j 


or more compactly. 



where 


c £ =(s £ )' e = dc £ e 5 -e T -dc E (d) r . 


( 21 ) 


( 22 ) 


(23) 


The relationships between field components given in the global, or beam coordinate system, and 
those in a system rotated by an angle 0 about the 3 direction (actuator system) are given by (see 
Figure 5) 


d = r £ d e = r £ e s = r,s t =(r/) t. 


where 



A 




X 


X 


d 2 

> E = < 

E 2 

► S = < 

s 2 

T = 

t 2 ► 


A. 


E 3 . 


X 




(similarly for D , E , S , and T ). 


(24) 


( 25 ) 
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The appropriate transformation matrices for this coordinate rotation are 


R = 


cos 2 0 
sin 2 0 

—2 cos0 sin0 


sin 2 0 
cos 2 0 
2cos0 sin0 


cos0 sin0 
-cos0sin0 > 
cos 2 0 - sin 2 0 




COS0 

-sin0 


0 


sin0 

COS0 

0 


0 

0 • 
1 


( 26 ) 


(27) 


In terms of the actuator coordinate system. Equation 22 then becomes 


D 


T 




(28) 


Substituting Equations 24 into Equation 28 yields constitutive relations expressed in terms of the 
global field variables. 


R/e s R £ R/eR s IfEl 
-R/e r R £ R/c £ rJ1s] 

(29) 

For convenience, electric fields and displacements will be defined as being specified along the 
actuator system 1 -direction only. The electric field within the piezoelectric material will also be 
assumed to be an average of the field strength between alternating electrodes. (See Reference 16 
for a complete description and discussion of the actual electric field distributions produced using 
IDE schemes.) Equation 29 is then 





T 


ef, eR s 11 eA 
-R/e r R/c £ R s JlsJ' 


(30) 


A material orientation angle that maximizes the actuator induced shear stress is desired here. This 
will occur for orientation angles of 0 = ±45° . For the case of 0 = +45° , the rotated material 
matrix of Equation 30 becomes 
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E n eR s 
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+ 2 c E + c E \ + c E 

T ^C ]2 T C 22 J T C 66 

— {c E + 2 c E + c E ] 
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4 V C 1 1 ^22 ) 
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4 V 1 1 22 / 


l(c E - 2 c E + C 
4 i ^ C i2 ~ 


Extracting the torsional terms yields 


A 

t 6 


1 

COI 
n to 

2(^33 e 3l) 

[£,1 

1 

1 

Ui 

1 

±(r £ -2c £ +c £ ) 

4 \ L 11 ^ c 12 T c 22 ) 

kl 


(31) 


(32) 


At this point , for simplicity, we will assume that the structurally effective portion of the blade is a 
closed, rectangular, thin-walled section (Figure 6). Following a mechanics of materials approach, 
the piezoelectric induced shear flow, q PE , , for such a structure may be written as 


Q PE 


= T t 

1 6 l PE 


(33) 


where t PE is the piezoelectric lamina thickness. The total piezoelectric torsional moment may 
then be expressed as 


Qpe 2hwq PE 2A e 2(^33 ^ 31)^1 ^pe' 

(34) 

where h and w are the height and width of the rectangular cross-section, and A e is the area 
enclosed by the wall centerline. 

Converting to standard piezoelectric coefficients, via Equation 33, yields 

Qpe = — (^33 ( C 11 — C 12 ) ^31 ( C 12 — C 22 ))^l *PE- 

(35) 

For further convenience, we will rewrite Equation 35 in terms of the assumed maximum 
produceable piezoelectric strain, A max , yielding 
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( 36 ) 
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max 


— ^33^1 


max 


Substituting Equation 36 into the expression for the generalized nondimensional piezoelectric 
control moments, from Equation 7, yields 


f 1 Qpe 

mil 2 R 2 


Jo 


1 A,c‘ 


mil 2 R 3 


(f 

1 _ 12 
7 E 

VV c n J 


(c t \ 


* 33 


V n 




' 22 




A max ^ PE 


OlfCc . 


(37) 


If the electrode layers are divided into P independently energizable sections, the average electric 
field intensity within the pth section is given by 


E lp ( V ,x)=£ w e p (\|/), 




p+1 


(38) 


where e p is the electric field generalized field strength, r p is the radial location of the inboard 

edge of the /?th electrode segment, and E { max is the maximum electrical field permissable without 
depolarization of the piezoelectric material. 


Substituting Equation 37 into 38 yields the generalized piezoelectric control moments in their final 
form 


f 1 Qpe 

Jo mil 2 R 2 


<& + m dx 


- 5 >J 
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(c E 

C 12 
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~ E 

C \\ ) 

“33 

L C 11 

C 11 J 


A ^t^dx. 


(39) 


Material properties appearing in Equation 39 for the piezoelectric laminae used in the numerical 
portion of this study are given in Table 4. 


3. NUMERICAL IMPLEMENTATION 

The equations of motion for the structural and aerodynamic degrees of freedom are rewritten in 
state variable form for numerical integration. The structural state derivatives thus will be of the 
following form: 
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struct 


-M'C 
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-M~'K 
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q 


struct 


+ 



fpE 

0 



with the structural state vector defined as 


q 


struct 



( 40 ) 


( 41 ) 


The definitions of M, C, K, etc., may be found in Appendix D. For each aerodynamic evaluation 
point along the blade, three state variables are required to define the lifting force and two to 
define the pitching moment. The aerodynamic state derivatives for each such point are 


/„r tJ + f 2i 
* 

f 2, 

* 

fu r 2i+/ 4i r 2l +/ 5l 

* 

r mil 

* 

[g u r m 2, + g 2; r 

m2 i + s 3 , 


with aerodynamic state vector 


q 


aero j 



r 2l . 



( 42 ) 


( 43 ) 


Expressions for the / and g coefficients appearing in the aerodynamic state space equations are 
given in Appendix E. 

The combined state vector is defined as 
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q = 


» aero 1 

q aero i\ J 


( 44 ) 


For N radial aerodynamic points, the resulting state space model will consist of a system of 
2 (L + M ) + 5 N first order differential equations with, in general, time varying coefficients. These 
state variable equations are numerically integrated using a fourth and fifth order Runge-Kutta- 
Fehlberg algorithm, with the integrations performed with respect to rotor azimuth angle, \j/. 

For the numerical case studies presented in this report, we have used one out-of-plane bending 
mode and three torsional modes. Here, (x) was defined as 


W { = x = x/R , 


which is the rigid body flapping deflection mode shape. 


(45) 


<D m (x) were assumed to be the first three torsional comparison functions developed by 
Karunamoorthy and Peters, 34 i.e., 


<J\ = 1, 

<I>3 = ^(30x 2 — 1 53c 4 — 7) 

(46) 

The odd numbered torsional comparison functions in Equation 46 have been found to be good 
approximations of the exact torsional nonrotating mode shapes for pinned blade-root boundary 
conditions. The even numbered polynomials approximate the cantilevered boundary condition 
nonrotating mode shapes. Use of polynomial approximations instead of the exact nonrotating 
mode shapes is done solely to simplify calculation of the integral coefficients appearing in 
Equations 6 and 7. 


4. RESULTS AND DISCUSSION 

4.1 Piezoelectric twist actuated rotor blade conceptual design 
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Three conceptual piezoelectric induced twist rotor blade designs were examined in this study. 
These designs were developed in order to illustrate the twist actuation capabilities of the three 
general cases of piezoelectric actuation suggested by inspection of Equation 36. These three 
general cases are: 

Case 1: d 3] * d 33 , c,f = c 22 • This is the case of actuation lamina possessing piezoelectric free- 

strain anisotropy, and in-plane stiffness isotropy. This corresponds to a configuration where 
actuation layers are composed of solid, or monolithic, PZT materials, and are polarized according 
to the IDE scheme. This case will be referred to as IDE/MON for the remainder of the 
discussion. 

Case 2: d M * d 33 , cf x * c 22 . This is the case where the actuation lamina possess both free-strain 

and stiffness anisotropy. This would be true of a piezoelectric fiber composite, interdigitated 
electrode actuation scheme. This case will be referred to as IDE/PFC. 

Case 3: d 3i = d 33 , c,f * c 2 E 2 . This is the case of free-strain isotropy (or near isotropy) but with 

stiffness anisotropy in the actuating layers. This would be the case for a piezoelectric fiber 
composite structure utilizing a conventional poling scheme, or a case similar to the Directionally 
Attached Piezoelectric scheme, originally proposed by Barret 12 and developed experimentally by 
Chen and Chopra. 13 This scheme will be refered to as DAP/PFC in the following sections. 

Stiffness, inertial, and actuation parameters for three conceptual helicopter blade designs, each 
representing one of the three cases above, were developed using full-scale helicopter parameters 
similar to those of the Sikorsky H-34 main rotor blade. The H-34 main rotor blade is a relatively 
simple, closed-section spar structure, and was easily idealized as a rectangular box section for the 
purposes of this study. These baseline blade parameters are summarized in Table 5. 

Structural parameters not identical to the baseline blade parameters are shown for each of the 
piezoelectric twist blade designs in Table 6. Uniform blade properties were assumed in each case 
for simplicity. The piezoelectric material thickness fractions, t PE „ given here were calculated 
assuming that the blade total mass of each design could be no greater than 120% of the baseline 
full-scale helicopter blade mass. (The choice of 120% was essentially arbitrary, but represents a 
reasonable weight constraint on the design of the conceptual piezoelectric twist blades.) As a 
result, the torsional natural frequencies of the blade structures vary somewhat from the baseline 
design. The resulting natural frequencies of all of the blade designs are given in Table 7. 
Aerodynamic parameters used for the numerical case studies were not varied between the designs 
and are shown in Table 8. 


4.2 Numerical twist actuation authority results 

Numerically generated twist actuation authority results for each of the three piezoelectric induced 
twist blade designs are shown in Figures 7-9. These results are for a typical one-g hovering flight 
condition, which corresponds to a thrust coefficient of Cj = 0.00465. One electrode segment 


19 



extending from x = 0,1 is assumed for all three structures. A sinusoidal electric field input with 
linearly increasing frequency and peak amplitude of E lmn was used to generate the frequency 

responses (amplitude and phase with respect to the electric field input signal) shown in the figures. 
In these figures the elastic twist is defined as the difference between the elastic torsional deflection 
at the blade tip and that at the blade root. 

Figure 7 illustrates the structural response of the IDE/MON case. A sustained oscillatory elastic 
twist magnitude of approximately ±1.25° to ±1.5° is generated for excitation frequencies below 
the first torsional frequency. At the first torsional resonance, which is predominately an elastic 
torsional response, the amplitude increases to approximately ±2.25°. A smaller torsional response 
occurs at the second and third torsional resonance frequencies. 

The resonant response at the second and third torsional frequencies was found to vary widely 
depending on the amount of material and aerodynamic damping present in the structure. As the 
torsional aerodynamic damping, from Equation 14, is in general proportional to b 2 , the 
corresponding aerodynamic damping for these two modes is almost negligible. Some additional 
form of damping is desirable then at these higher frequencies to avoid unrealistically large 
torsional deflections. As such, a level of material damping equivalent to 0.5% of critical damping 
was assumed for each of the cases presented here. 

The actuation results for the IDE/PFC lamina design are shown for the same flight condition in 
Figure 8. A level of actuation capability on the order of ±1° to ±1.25° of elastic twist below the 
first torsional resonance, and ±1.5° at the first torsional mode resonance frequency is shown here. 
This is a level of performance slightly less than that demonstrated with the IDE/MON 
configuration. Although this may seem to imply that monolithic PZT laminae are more desirable 
for inclusion in piezoelectric actuated structures, manufacturing and poling nonplanar composite 
structures with solid PZT layers may not be practical. Piezoelectric fibers on the other hand could 
be incorporated into complex composite aerospace structures using, for the most part, established 
fiber composite construction techniques. 

Figure 9 displays the twist actuation capabilities of the DAP/PFC blade design. Structurally, the 
DAP/PFC blade is identical to the IDE/PFC blade design, although the DAP/PFC blade utilizes 
conventional poling of the piezoelectric fibers. Relatively low nonresonant twist actuation is 
demonstrated for this actuation case, i.e., around ±0.2° to ±0.25° of elastic twist. 

Comparison of the elastic twist actuation response of all three cases is shown in Figure 10. The 
effect of the large free-strain anisotropies present in the IDE schemes on the magnitude of elastic 
twist is readily apparent. Both IDE poling cases exibit generally four to five times the twist 
actuation magnitudes of the conventionally poled configuration. Such magnitudes of elastic twist 
are generally regarded as being sufficient for practical use in a vibration reduction scheme using 
individually controllable blade twist. 
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5. CONCLUSIONS 


A simple helicopter rotor blade aeroelasticity analysis was developed and used to numerically 
demonstrate the twist actuation potential of embedded piezoelectric actuators for three nominally 
full-scale helicopter rotor blade designs. It was numerically demonstrated that useful nonresonant 
levels of oscillatory blade twist, i.e., on the order of ±1°, can potentially be produced without the 
addition of an excessive amount of piezoelectric actuator mass or saturation of the piezoelectric 
actuator materials, using an interdigitated electrode poling scheme with either a piezoelectric fiber 
composite or monolithic PZT actuation design. 

The analysis and numerical model in its present form (i.e., with rigid flapping, elastic torsion and 
stall aerodynamics), should be sufficient for an examination of the potential of piezoelectric twist 
actuation to alleviate high oscillatory control loads induced by blade stall flutter. 35 Such a study is 
underway by the authors. Improvements to this model, such as the addition of multiple flapwise 
bending modes and a simple nonuniform inflow model, are also being undertaken. 
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APPENDIX A. COORDINATE SYSTEMS AND TRANSFORMATIONS 


Several coordinate systems are used to describe the flap-torsion kinematics of the rotor 
blade. These coordinate systems are shown in Figure A1 , and described below. 

1. XjYiZi is an inertial reference frame with origin at the hub center of rotation, Xi axis 
oriented in the direction of forward flight, and negative Zi axis aligned with the rotor 
rotation vector. 


2. X]Y;Zi is a hub-fixed rotating reference frame, with the Zi axis in the direction of the 
rotor rotation vector. The X]YiZi system is rotated about the Zo axis by the blade 
azimuth angle \\f. The Xi axis is coincident with the undeformed blade elastic axis. 

3. x 0 yoZo is a blade-fixed reference system located at an arbitrary point (x) along the blade 
undeformed elastic axis. 

4. xyz is a reference system fixed in the deformed blade, and translated with respect to 
the x 0 yoZo system. 

5. x 2 y 2 z 2 is an intermediate system obtained by rotation of the xyz system about y by the 
angle w' . 

6. X 3 y 3 Z 3 is a reference system fixed in the deformed blade obtained by rotating the x 2 y 2 z 2 
system about y 2 by the total rotation angle 0—0 cgn + tf) . x 3 is tangent to the deformed 

elastic axis. 

Coordinate transformations due to rotations between the various systems are given below. 

Transformation from the inertial system (X]YiZi) to hub-fixed rotating system(Xi Y {Z\): 

-cost)/ sinvj/ 0 

'T 0 = siny cost]/ 0 

° 0-1 

(Al) 


Transformation from the blade-fixed xyz system to the intermediate deformed system 
x 2 y 2 z 2 (rotation about y by w ' ): 


cosw' 0 sinw' 

2 T,= 0 1 0 

-sinv/ 0 cosw' 

(A2) 
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Transformation from the x 2 y 2 Z 2 system to the x 3 y 3 z 3 system (rotation about y 2 by 0): 


1 0 

3 T 2 = 0 cos0 
0 -sin0 


0 

sin0 

COS0 


(A3) 
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APPENDIX B. GENERALIZED FORCE DISPLACEMENT RELATIONS FOR A 
COMPOSITE THIN WALLED BEAM STRUCTURE 


Rehfield 36 presents a general theory for characterizing the structural properties of thin 
walled composite beam structures. Portions of this theory necessary for deriving the beam 
stiffness properties of an idealized rectangular thin walled rotor blade structure are given 
below. 


Membrane resultant stresses and strains for a composite laminate are related by the 
membrane stiffness matrix, A, as follows, 




1 

> 

> 

> 

Os 

I 


X ' 

n 2 

> - 

A 2 ^22 ^26 


s ° 2 > 

" e . 


_^16 ^26 ^66 _ 




(Bl) 


The elements of A are formed from the individual laminate plane stress stiffnesses 

rr^l 37 • 

according to Classical Laminated Plate Theory , i.e., 


=XL £ ). 


v k ’ 


(ij = 1,2,6), 


k = 1 


(B2) 


where n is the number of plies of the laminate. 

With the assumption that the hoop stress, N 2 , is zero, the hoop strain, S 2 , may be written 
as 


S 


0 

2 



(B3) 


Applying equation B3 to B 1 , and rewritting the extensional and shear stresses in terms of 
the remaining strains yields 



(B4) 


where the K stiffnesses are given by 
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A 2 ^26 
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(B5) 


(B6) 


(B7) 


The relation between the generalized elastic beam forces and generalized displacements 
used in the blade equations of motion of this paper may then be expresed as 
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C» 

Cy 4 

C 15 

M x 
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Q 4 

C 44 

C 45 



c a 

C 45 

c 55 


(B8) 


The elements of matrix C correspond to the beam stiffnesses of the blade equations of 
motion as follows: 


C u = ^ K n ds ~ EA , 

(B9) 

4 A 2 f 


C M = — j—\K 22 ds == GJ , 


c J 

(BIO) 


Q; A, z 2 ds~EI m , 



(Bll) 

2 A r 


C„= 'U n ds. 


c J 

(B12) 


Integrations in equations B9-12 are performed around the contour of the beam section, 
with A e the area enclosed by the contour perimeter, and c the contour length. 

For the blade structural configurations discussed in this paper, coupling stiffnesses C 15 and 
C 45 are identically zero. The extension-torsion coupling stiffness, C 14 , is zero for the 
monolithic PZT actuation structure. The fiber composite structures will inherently contain 
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some extension twist coupling, resulting in a nonzero C14. This coupling, however, is 
relatively small for the particular configurations studied in this paper and was not included 
in the blade equations of motion. Future, more extensive, developments of the blade 
equations will include these elastic couplings. 
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APPENDIX C. DEFINITION OF SECTION ANGLE OF ATTACK AND 
SECTION ROTATION RATE WITH RESPECT TO THE AIR MASS 


C.l Section angle of attack: 


The two dimensional angle of attack of the ith airfoil section (a,) is defined by 

-if U A 

a = tan — - 

(Cl.l) 

where U x and Uy (omitting subscripts) are the components of the nondimensional section 
airflow velocity normal and parallel to the section chord. U x and U Y expressed in the 
customary airfoil coordinate system (Figure Cl) are 


U x 




with the total nondimensional section airflow velocity given by 

u = Ju^ r +u/ . 


(Cl. 2) 


(Cl. 3) 


Spanwise flow, i.e., flow in the x 3 direction is neglected. 

V z and V y are the components of the section airfoil velocity with respect to the air mass in 
the x 3 y 3 z 3 system (see Appendix A), viz. 


\ = <V y 
V. 


= V - 3 V„ 


= 3 T 2 2 T, I T 0 °V a - 3 T 2 2 T 1 1 


(Cl. 4) 


°V a and *V b are, respectively, the velocity of the air mass in the inertial system, and the 
velocity of the airfoil section due to blade motion with respect to the XoyoZo system, given 
by 
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-V x cosa / 


0 

II 

a 

> 

o 

0 

ii 

rQ - 


V 

k J 


rw' 


(Cl. 5) 


Making appropriate small angle assumptions for a s , w ' , and 0 (and implicitly for 0 con and 
<|)), yields for V: 


X 


V m cost)/ - vw' 


jcw'vv' 

Vy 

> = •< 

—V" sin V(/ - V^w'Q cosxj/ - V0 

> — « 

xQ + xvv'0 > 



sin\)/ - V^w'cosx)/ -v 


— xQ0 + xw' j 


V m cost)/ — Vvv'jc — xw'w' 

—V„ (w'0 cosx|/ + sinxj/)-v0 - x(Q + vv'0) > 
^(©sinx)/ - w'cosxjr)-v + jc(Q 0 - vv') 


(Cl. 6) 


Dividing by OR yields the components of the nondimensional section velocity in the blade 
coordinate system (3-system) 


U, 

V> 

u. 


pcosxj/- 'k TPP w + +w + w + x 


-jlsinxj/ - |iw + 0cosx|/ -’k TPP Q - x 


( 


1 + w + 0 


V 


p.0 sinxy - pw + cost]/ -X TPP + x 0 - w + 


(Cl. 7) 


where (*) denotes differentiation with respect to the rotor azimuth, and ( + ) differentiation 
with respect to x/R. Substituting U z and U y into equation Cl.l and then into equation 
Cl. 2 yields the section angles of attack. 

C.2: Section rotation rate: 


The rotation rate of the airfoil section with respect to the air mass will be defined as 

* d F * * 

E = = £e +£*', 

d\\f 

(C2.1) 
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♦ 4 

where Ee is the rotation rate of the airfoil due to blade pitch motion, ande w - is the 
effective rate rotation of the airfoil due to the component of the rotor rotation vector lying 
along the section pitch axis, x 3 . 


Ee is given simply by 


Be 


30 * * * 

r— =e = e«*+4>. 
3\J/ 


(C2.2) 


Ew' may be found by examination of the vector component of the rotor rotation lying 
along the x 3 , i.e.. 


£ 


w" 


dt 3\|/ 


= Q£w' = £2k, i 


3 » 


(C2.3) 


where k| is a unit vector in the direction of the positive Zi axis, and i 3 a unit vector in the 
direction of the positive x 3 axis. Rewriting C2.3 yields 


Eh-- = k, i, = sin w' = w + , 


for small rotations. 


(C2.4) 


The total rotation rate with respect to the air mass may then be written as 

* * * 

E = 0 con + <|)+ W + . 


(C2.5) 
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APPENDIX D. SYSTEM MATRICES 


Equations 6 and 7 may be rewritten as a system of M + N ordinary differential equations 
of the form 


M x+ Cx+ Kx = f„„, + f„ + f, , 


(Dl) 


where x is a vector of structural generalized displacements, given by 
x = {vv, vv 2 ••• w L <j), <{) 2 ••• <|> m } T . 

(D2) 


The symmetric mass matrix, M, in equation Dl is given by 


M = 


[/°°] [J 00 ] 

[j°°] T [/i: 00 ] ’ 


with 



(D3) 


(D4) 

(D5) 

(D6) 

(D7) 
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Tf (X) 
A 11 

Tf 00 
A 12 

Tf 00 

A 1A/ 

II 

8 

* 

rr 00 
A 21 

Tf 00 
A 22 

* 


Tf 00 

_ A M1 



(X) 

K MM 

*=-£ 

* o 

gl m n 

dx. 


C is a symmetric structural damping matrix. 


with 


C = 


M [o] 

[0] [•*"] 


l n = 


22 

11 

22 

21 


22 

12 

22 

22 


f 22 

1 L\ 


1 22 
l \L 


* x22 

1 LL 


*/ 22 = f 1 ^ ^pi 

mn J 0 


1 milR 


A W m ++ W„cE , 

i4 m n 7 


= 


* rri 1 

mn 


* wf 11 

A 11 

’AT, 

* If 1 1 
A 21 

'K, 

* T^ll 
A W1 

... 

= f‘_ 

gV 

J ° mQ/? 4 


* jt^ll 
A 1A/ 




MM . 


®_ + 0 + <£c. 


K is a linear stiffness matrix of the form 


(D8) 

(D9) 


(DIO) 


(Dll) 

(D12) 


(D13) 

(D14) 
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where 


(D15) 


K = 


[/-] ir°] 

[,/ ,0 ] r [K**] 


jWW _ 


rf r /ii+ o / , i 

^ 'ii — i n 

+ /?) 

... 

/r/u o/ii 
V ' 1L ' 'li. 

+ «)' 

( t i u +°i u 

^ *L1~ 1 LI 

+ tf) 

... 

f r /"+°/ n 

V ' LL T 1 LL 

+/2). 


T I l m [ = f 1 — T—w+w+dx. 

mn J °mQ 2 /? 2 


0 rll _ 
1 mn 


K a 


Q 2 R 2 f mdx 
Jo 

i EL 


^ + (0)W„ + (0), 


7 22 ■ f "7 1 4 W m ++ W n ++ dx , 
mn J °mQ 2 i? 4 m " 


y 10 = 
1 J 


10 


i , 'ii 


1 J L 1 


y 10 

1 J \M 


J'° 

1 J LM 


(D16) 

(D17) 

(DIB) 

(D19) 


K™ = 


,JZ s few +<!><£?, 

l mn Jq m n 7 

/ 1^00 I Tr^OO . j /- 1 1 , 0 rrlX) \ /e^OO.T’k^OO, 1 , 0 r^OO \ 

^Aii + A„ +Au+ A u J ••• [R lM + + A 1M + K lM ) 

/^00.7T^00.jrll.0K^00\ / jv" 00 , T 00 , jr^ 1 1 , 0^00 \ 

a:“ = f 1 —— r- o o dj , 
mw « 0 ^ 2 m n 

t k u s f 1 5 

mn Jo D 2 Z > 2 


(D20) 

(D21) 


f i r k 2 + _ + _ 

, , ^VO m + O n dJ, 

'mQ 2 /? 2 fl 2 


(D22) 

(D23) 

(D24) 
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(D25) 




f Mro , f PE and / 0o are, repectively, the nondimensional aerodynamic forces, 

nondimensional piezoelectric induced twisting moments, and nondimensional steady state 
inertial propeller moments, and are shown below: 



(D27) 

(D28) 

(D29) 



(D30) 

(D31) 
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(D32) 

J° E 

m 

e |W m <£t , 
Jo m 

(D33) 

,7i = 

\'exW m + dx, 
Jo m 

(D34) 

a: 0 = 

fi k 2 


I " & m dx. 


m 

lo R 2 

(D35) 
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APPENDIX E. AERODYNAMIC STATE SPACE COEFFICIENTS 


Coefficients associated with lift formulation: 



i 

ii 

(El) 


f 2 =(aX 1 U y +be'S z )/k 

(E2) 


A—M,wjk 

(E3) 



(E4) 

fs=f 4 

|^(/ X AC Z + AC Z + ^ U y J 

(E5) 

Coefficients associated with pitching moment formulation: 



=~ a J k 

(E6) 


s 2 =-r m / k2 

(E7) 


gj = g 2 UAC m -E m U Y /k. 

(E8) 


38 



Table 1 . Ordering of nondimensional parameters appearing in equations of motion. 

(e = 0.10). 


u/R = 0( e 2 ) 

El/mQ 2 R 4 =0(e 3 ) 

vjR = 0(e 2 ) 

G/mfi 2 = 0(e -6 ) 

w/R = 0(e) 

EA/mQ 2 R 2 = 0(e“ 2 ) 

<t> = 0(e) 

t'/R 2 = 0(e 4 ) 

0 = 0(e) 

^/* 2 = 4 = 4 ) 

Hi 

II 

"ST' 

u> 

5;/A^ 4 =0(e 4 ) 

II 

'TrT 

£ 2 ‘/AR 3 =0(e 3 ) 

x/R = 0(1) 

y//? 4 =0(e 9 ) 
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Table 2. Lift based ONERA stall coefficients for the OA212 airfoil. 


Parameter 

Value 


Physical Description 

X 

0.2 


time delay parameter 

s 

5 

5 

fa - 4 < 1 + U3AC <) 


apparent mass quantity 
lift coefficient to pitch rate 
relation 

W 

0.10+ 0.023(|a| - 13°)w(|a| - 13°), 
0.3, 

|a| < 21.7° 
|a| > 21.7° 

damping factor 

(m is a unit step function.) 

d 

0.105/w 


stall natural frequency 

e 

2 -5.1tan“ 1 [l.2l(|a|-13°)]H(|a 

|-13°) 

phase shift parameter 
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Table 3. Coefficients for stalled pitching moment circulation equation. Values are for the 

"mean airfoil" of Reference 27. 


Parameter 

Value 

To 

0.2 

r 2 

0.2 

ao 

0.25 

a 2 

0.1 

e 2 

0.573 
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Table 4. Lamina structural properties. Piezoelectric lamina properties are derived from 

values reported by Rogers and Hagood. 


Property 

passive blade 
structure 

monolithic PZT 
lamina 

IDE/PFC 

lamina 

DAP/PFC 

lamina 

Cn 

12.8e6 psi 

9.6e6 psi 

4.5e6 psi 

4.5e6 psi 

C22/C1I 

1 

1 

0.60 

0.60 

C 12 /C 11 

0.375 

0.29 

0.24 

0.24 

C66 

4.0e6 psi 

3.4e6 psi 

8.3e5 

8.3e5 

d 3 i/d 3 3 

- 

-0.5 

-0.4 

1 

p 

0.1 lbs/in 3 

0.27 lbs/in 3 

0.21 lbs/in 3 

0.21 lbs/in 3 

max 

- 

500 pe 

500 p£ 

250 jxe 
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Table 5. Baseline helicopter rotor parameters. 


Parameter 

Baseline blade value 

12 (rad/sec) 

23.2 

R (in) 

336 

y 

9.44 

c/R 

0.0488 

n 

4 

a 

0.0622 

Cj 

0.00465 

f/nR 2 

0.015 

m e = Vk 9 /Q 2 / 0 

16.0 

m (slinch/in) 

0.00142 

m/c 2 (slinch-in 2 /in) 

0.0175 

U/h 

0.000327 

e/c 

0 

GJ 

0.00552 

I p Q 2 /? 


h (in) 

1.84 

w (in) 

6.11 

t (in) 

0.141 
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Table 6. Structural parameters for numerical examples. 


Parameter 

IDE/MON 

IDE/PFC 

DAP/PFC 

GJ 

7 p Q 2 /? 

0.00447 

0.00365 

0.00365 

Y 

8.28 

8.28 

8.28 

tpE A 

0.1875 

0.300 

0.300 
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Table 7. Calculated nondimensional blade natural frequencies for numerical examples. 


Parameter 

Baseline blade value 

IDE/MON 

IDE/PFC 

DAP/PFC 

(b p 

1.0 

1.0 

1.0 

1.0 


6.10 

5.60 

5.14 

5.14 


18.19 

16.64 

15.22 

15.22 

03 

61.02 

59.67 

58.58 

58.58 
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Table 8. Aerodynamic parameters for numerical examples. 


Parameter 

value 


Cl 

-|(1+1.4M’) 


C 2 



C3 

'Ijr 

-—(-1.26 — 1.53 tan ■' (1 5 (M - 0.7))) 

M 

0.30 


y ac 

0 


N 

5 


X 

{0.28 0.44 0.60 0.76 

0.92} 

A 

{0.16 0.16 0.16 0.16 

0.16} 
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Figure 1. Conceptual diagram for a unified rotor blade piezoelectric aeroelastic analysis. 
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* 


3 



layers, top and bottom 


Figure 4. Material axis system for piezoelectric fiber composite assuming interdigitated 

electrode poling scheme. 
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Figure 5. Rotation of piezoelectric fiber composite material system with respect to the 

global (blade) system. 
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blade cross section: 



structurally effective cross section: 

z 3 



Figure 6. Idealized rectangular, thin-walled, closed-section piezoelectric blade structure. 
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magnitude, deg 




Figure 7. Blade frequency response: monolithic PZT (MON) twist actuation with 
interdigitated electrode poling scheme; C T = 0.00465, p = 0.0. 
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magnitude, deg 




Figure 8. Blade frequency response: PFC twist actuation with interdigitated electrode 

poling scheme; C T = 0.00465, p = 0.0. 
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magnitude, deg 



frequency, P 



Figure 10. Comparison of elastic twist frequency response for IDE/MON, IDE/PFC, and 
DAP/PFC actuation schemes; C T = 0.00465, fi = 0.0. 
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Figure A1 . Rotor blade coordinate systems. Note that the blade undeformed elastic axis 
lies along Xi. (Section pitch angle, 0, is not shown.) 


57 





REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


Prtiic reporting burden for (hit colectior of Information fc estimated to average 1 hour per reeponao, including the lime lor re viewin g instructions, searching relating date sources, 
gathering and maintaining the data needed, and completing and reviewing the cotiection of In fo r ma tion. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headq ua rters Services, Directorate for infornwtion Operation s and Reports, 1215 Jefferson Davis 
Highway, Suite 1204, Arlington. VA 22202-4302. and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188). Washington. DC 20603, 


1. AGENCY USE ONLY ( Leave blank) Z REPORT DATE 


May 1996 


4. TITLE AND SUBTITLE 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


An Aeroelastic Analysis of Helicopter Rotor Blades Incorporating 
Piezoelectric Fiber Composite Twist Actuation 


505-63-36-02 


6. AUTHOR(S) 

W. Keats Wilkie 
K. C. Park 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRES$(ES) 

NASA Langley Research Center 

Hampton, VA 23681-0001 
Vehicle Structures Directorate 
U.S. Army Research Laboratory 

NASA Langley Research Center, Hampton, VA 23681 -0001 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 and 
U.S. Army Research Laboratory 
Adelphi, MD 20783-1145 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


10. SPONSORING / MONITORING 
AGENCY REPORT NUMBER 


NASA TM-1 10252 

ARL-MR-328 



12a. DISTRIBUTION /AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 05 


13. ABSTRACT (Maximum 200 words) 



A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, 
interdigitated electrode blade twist actuators is described. The analysis consist of a linear torsion and flapwise 
bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin 
procedure is performed upon the rotor blade partial differential equations of motion to develop a system of 
ordinary differential equations suitable for numerical integration. The twist actuation responses for three 
conceptual full-scale blade designs with realistic constraints on blade mass are numerically evaluated using the 
analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to 
two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling 
configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of 
a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated 
electrode actuators was found to be four to five times larger than that obtained with the conventionally poled 
actuators. 


14. SUBJECT TERMS 


Helicopters, Piezoelectric Actuation, Aeroelasticity, Adaptive Structures 


15. NUMBER OF PAGES 


116. PRICE CODE 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rav. 2-89) 
Prescribed by ANSI Sid. Z38-18 
298-102 
















Wilkie, W. Keats, and Park, K. C., " An Aeroelastic Analysis of Helicopter Rotor Blades Incorporating 
Piezoelectric Fiber Composite Twist Actuation," NASA TM 110252, U.S. Army Research Laboratory Report 
328, May 1996. 

ERRATA (compiled as of 7/23/96) 

p.5: (clarification) Although not stated, small angles are assumed implicitly for 0 CO „ and 4>. 

p.5: In Eq. (1) several subscripts were omitted, and some exponents arc in error. 




+ 


El 


\mil 2 R 4 


vc 



£7 m 

jr 

milR 



K 

mil 2 R 


p.5: Eq. (2) has some missing subscripts and erroneous exponents. The equation should read 

-ML-,-)' V _^L__r_as_V hLe 

) {milR ' 9 ) mflV R 2 ~ 

p.6: Several exponents and subscripts in Eq. (6) arc in error. The equation should read 

iiX'j'wWJx 
+j\x, \'!L\wrw;'dx 

ft Jo mClR' ' " 

l ^ T i El K 

= f-^jzrW.dx-CO'.. f 'cW.dx + 9 m feTW/dt) 

J °mCl 2 R \ h ’ m Jo * ) 

p.7: Eq. (7) has missing subscripts and erroneous exponents. The equation should read 


m~ 1 


+ y i> m 


1 mtlR 
r'k. 2 


Jo 

+ Jf, ]W,<I>,dF+ Yw, jWl^de 

= f -pTZi9 , d x + f'- 

] Omil 2 R 2 r M/nfl’S' ' ( ""/Jo fl 2 ' 


p.9: (clarification) Eq. (13) is an approximation valid for f/ K < < U x , and small 0. The exact expression is: 

£Z{[ f; r,H 0 ». + *)(,-, -f/ y .sin(^+^)| f ](r, +r j ) + /T Oi cos(0_+^)| f }JJ" W^dr 


1 



p.10: Subscript on the fourth term of line two is missing. Change to read 






p.15: (clarification) In Eq. (33), T 6 should be the piezoelectric actuation shear stress, i.e., the stress due 

solely to an application of an electric Held with no applied mechanical strains. Also, r in Eq. (34) is the 
perpendicular distance from the tangential direction on the section contour to the origin. 

p.15: The equation reference immediately above Eq. (35) is in error. The correct reference is to Eq. (23). 

p.16: Exponents in the first terms of Eq. (37) and Eq. (39) arc in error. Both terms should be 

" r '-Qs-r*:* ..." 

p.17: In Eq. (40), a subscript is incorrect in the second term after the equals sign. This term should read 

"... + ?/*+?._} . " 

p.33: A subscript is missing on the last term of Eq. (Dl). The equation should read 

M x+ C x+ Kx - f -iW + 1 pe + fff 

p.35: The exponent on R in the denominator of Eq. (D18) is incorrect. The equation should read 

= — 771 ^*(0R‘(0) 

Q 2 «’|md5r 

Jo 

p.36: The exponent on R in the denominator of Eq. (D26) is incorrect. The equation should read 

‘’O— 77T *„(0)*.(0) 


Jo 


mdx 


p.36: The sentence following Eq. (D26) should be changed to read: 

* f arro , fp E an d f 9 are, respectively, the nondimensional aerodynamic forces, nondimcnsional 
piezoelectric moments, and nondimcnsional blade control pitch moments, and are shown below:” 

p.36: Incorrect exponent in Eq. (D31). The equation should read 

Q - 9s. 


p.37: Missing subscripts in Eq. (D32). The equation should read 



jy 



I 

II 

J°L 


A 

K 


K 


K 


2 



